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1.0  INTRODUCTION 


This  report  contains  a  theory  of  the  anomalous  second  harmonic  generation  (SHG)  ob¬ 
served  in  certain  atomic  vapors  [Refs.  1-9]  and  a  theory  of  the  band  structure  for  su¬ 
perlattices/  multiple -quantum-wells.  The  first  theory  suggests  that  it  may  be  possible 
to  double  the  1.3 n  radiation  at  high  intensities  and  large  outputs,  without  considera¬ 
tion  of  damage  inherent  to  nonlinear  crystals.  The  second  is  a  preliminary  background 
for  investigations  into  radiative  interactions  within  nonlinear  heterostructure  crystals 
which  have  several  orders-of-magnitude-larger  nonlinear  susceptibilities  than  homoge¬ 
neous  crystals. 

The  theory  of  the  anomalous  SHG  is  based  on  global  symmetry  breakdown  in  an  ini¬ 
tially  isotropic  gas  resulting  from  the  off-diagonal  order  created  by  the  fundamental  co¬ 
herent  mode  in  collision-induced  radiative  coupling  channels.  The  effect  is,  therefore, 
intimately  connected  with  the  modification  of  collision  dynamics  under  applied  laser 
fields.  In  the  new  theory,  collisional  and  radiative  processes  are  taken  into  account  si¬ 
multaneously.  It  is  based  upon  the  Hilbert  space  of  the  representation  as  it  is  defined  by 
the  scattering  states  of  pairs  of  colliding  atoms.  It  is  an  alternative  to  the  phenomeno¬ 
logical  theories  used  h  pressure- induced  resonances  in  four- wave- mixing  [Refs.  10  -  11]; 
to  the  semiclassical  theories  used,  in  collision-induced  absorption  [Ref.  12]  and  true  and 
avoided  crossings  in  intense  fields  [Ref.  13],  and  to  the  projection  algebra  of  the  density 
matrices  used,  for  example,  in  the  collisional  redistribution  of  light  [Ref.  14]. 

In  Section  2.1  the  experimental  background  of  anomalous  SHG  is  described.  In  Section 
2.2  the  dynamics  of  a  single  collision  complex  in  a  coherent  radiation  field  are  investigated 
using  the  new  theoretical  approach.  This  section  shows  that  the  Mollow  spectrum  as 
well  as  the  spectrum  of  the  Rabi  oscillations  is  modified  by  collisions.  In  Section  2.3 
the  problem  of  N  collision  complexes  and  the  quantum  description  of  phase-matching  is 
considered.  In  Section  2.4  an  expression  for  the  macroscopic  second  harmonic  polariza¬ 
tion  vector  is  derived.  In  Section  2.5  a  specific  gas  mixture  which  might  be  useful  for  the 
doubling  of  the  oxygen-iodine  frequency  is  discussed.  In  Section  3.0  the  band  theory  of 
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superlattices/multiple-quantum-wells  is  discussed.  It  is  formulated  in  the  Bloch  repre¬ 
sentation  of  one  of  the  constituent  materials,  and  spontaneously  exhibits  the  symmetry 
breaking  along  the  superlattice  axis,  which  causes  the  folding  of  the  original  Brillouin 
Zone  onto  a  thin  region  around  the  center  of  the  zone.  A  single  homogeneous  crystal 
band  splits  into  many  subbands.  The  series  expressions  for  the  energies  of  these  sub¬ 
bands  is  given  and  the  radiative  selection  rules  are  discussed. 
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2.0  SECOND  HARMONIC  GENERATION  FROM  GASES 


2.1  EXPERIMENTAL  BACKGROUND 

The  observation  of  SHG  from  an  isotropic  atomic  vapor  was  first  reported  in  1978  (Ref.  V. 
In  the  experiment  atomic  T1  vapor  was  subjected  to  external  dc  fields.  However,  there 
was  SHG  from  the  T1  vapor  even  when  all  the  external  fields  were  turned  off,  and  when 
u  was  tuned  to  half  the  transition  frequency  between  the  6s  Pi.  and  7 2  Pi  levels.  The 
experimenters  had  used  an  JV2-  pulsed  dye  laser  which  produced  7  nsec  pulses  at  a 
fundamental  wavelength  Aw  =  5860  A  (with  FWHM  1  to  10  GHz).  The  T1  vapor  was 
in  a  heat  pipe  at  a  temperature  of  860  0  C.  The  second  harmonic  signal  was  at  A2u  = 

2930  A,  and  its  polarization  followed  the  polarization  of  the  fundamental  mode.  The 
second  harmonic  intensity  /2u  scaled  as  N3  ,  where  N  is  the  density  of  the  vapor  and 
Iu  is  the  intensity  of  the  fundamental  mode.  There  was  a  narrow  frequency  interval  in 
which  SHG  occurred.  The  SHG  was  destroyed  by  small  detunings  of  about  2  GHz  on 
either  side  of  the  resonance.  When  an  external  magnetic  field  was  applied  to  the  va¬ 
por,  it  had  no  effect  on  the  second  harmonic  signal  A2w  =  2930  A.  This  ruled  out  Ml- 
type  transitions.  The  experimenters  advanced  the  hypothesis  that  multiphoton  ion¬ 
ization  of  T1  atoms  and  a  subsequent  charge  separation  between  faster  electrons  and 
slower  ions  might  be  responsible  for  the  anomalous  SHG.  If  E0  was  a  nearly  static  field 
arising  from  multiphoton  ionization,  then  it  could  cause  SHG  via  four-wave-mixing: 

P2uL  =  X(3)  :  E„EUE0.  However,  this  hypothesis  was  not  very  satisfactory,  since  it 
implied  that  the  scaling  relationship  between  I7u  and  Jw  is  I2u  oc  Ilk+2,  where  k  is  the 
minimum  number  of  photons  required  for  the  ionization  of  an  atom  by  the  fundamental 
mode.  For  T1  atoms,  k  >  3.  Thus,  the  multiphoton  ionization  scaling  was  in  contradic¬ 
tion  to  the  observed  scaling. 

The  ionization  model  was  tested  in  a  later  experiment  [Ref.  3].  The  atomic  Na  vapor  in 
a  heat  pipe  was  irradiated  with  intense  Aw  =  1.064^m  pulses  of  28  psec  duration  and 
a  weak  SHG  at  A2w  =  5320A  was  observed  .  The  wavelength  Aw  =  1.064/i  is  not  reso¬ 
nant  with  any  of  the  Na  transitions.  But  the  experimenters  estimated  that  for  intensities 
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above  5  x  101  lW j cm7,  the  fundamental  beam  would  cause  partial  ionization  of  the  va¬ 
por  due  to  multiphoton  transitions.  A  field-gradient  mechanism  was  suggested  to  explain 
the  SHG  from  the  original  isotropic  vapor.  In  the  model,  the  influence  of  the  multipho¬ 
ton  ionization  was  secondary;  it  augmented  the  SHG  caused  by  the  field-gradient  mech¬ 
anism.  The  latter  implied  that  the  dc  electric  fields  in  the  gas  were  on  the  order  of  50 
V /cm,  and  that  the  second  harmonic  intensity  scaled  as  I* . 

The  two  experiments  yielded  significantly  different  results.  The  second  harmonic  inten¬ 
sities  observed  under  picosecond  pulses  in  Na  were  several  orders  of  magnitude  smaller 
than  those  in  Tl.  The  I3u  had  completely  different  scaling  properties  with  respect  to  the 
fundamental  beam  intensity  in  the  two  experiments.  Such  facts  raised  strong  doubts  that 
the  SHG  had  the  same  underlying  mechanism  in  the  two  experiments.  These  questions 
were  investigated  in  another  set  of  experiments  [Ref.  6].  The  experimenters  used 
atomic  Na  vapor  in  a  cell  and  a  pulsed  dye  laser  whose  output  could  be  tuned  to  some 
of  the  Na  transitions.  The  5P  -  3S  and  4D  -  3S  transitions  yielded  SHG,  but  the  5S  - 
3S  transition  did  not  yield  any  SHG.  The  SHG  was  tunable  about  the  resonances.  The 
shapes  of  the  tuning  curves  were  sensitive  functions  of  the  vapor  density.  This  suggested 
that  the  SHG  was  caused  by  a  bulk  interaction.  Simultaneously  with  the  SHG,  there 
was  a  broad-band  incoherent  fluorescence  at  2 w.  To  test  for  ionization  effects,  the  ex¬ 
perimenters  irradiated  the  cell  with  intense  pulses  from  a  Nd-YAG  laser.  The  pulses 
had  no  effect  on  the  SHG.  This  ruled  out  multiphoton  ionization  as  the  mechanism  for 
the  anomalous  SHG.  The  experimenters  introduced  a  higher  mode  structure  into  the 
fundamental  beam  in  order  to  test  for  field-gradient  effects.  This  also  had  no  effect  on 
the  SHG  and  ruled  out  the  field-gradient  mechanism.  Nor  could  their  results  be  ex¬ 
plained  by  quadrupole  transitions,  because  these  yield  quantum  efficiencies  that  are  six 
to  eight  orders  of  magnitude  less  than  the  experimental  efficiencies  [Refs.  2,6].  The  ex¬ 
perimenters  varied  /„  from  105 W/cm2  to  10®  W/cm2.  For  the  S  -  D  transitions,  they 
found  that  I3v  scaled  as  at  low  intensities.  But  I2u  rapidly  began  to  scale  as  1 2,  and 
then  as  Iu ,  as  Jw  was  increased  to  the  values  in  the  upper  part  of  its  range. 

The  anomalous  SHG  has  been  observed  in  other  atomic  vapors:  Li  [Ref.  7],  Ba  [Ref.  4], 
K  [Ref.  9],  Hg  [Ref.  5],  and  Zn  [Ref.  8].  The  experiment  with  Zn  is  of  special  interest  in 
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that  an  attempt  was  made  to  measure  the  effective  electric  fields  which  could  generate 
the  observed  second  harmonic  fields  obtained  from  the  same  atomic  levels.  The  electric 
fields  were  5  x  103  V/cm  and  4  x  103  V/ cm  for  the  D  -  S  and  S  -  S  transitions,  respec¬ 
tively.  These  values  are  two  orders  of  magnitude  larger  than  the  value  inferred  in  Ref. 

3. 

2.2  A  NEW  FORMULATION  OF  THE  RADIATIVE  COUPLING  TO  A  COLLISION 
COMPLEX 

This  section  presents  the  results  of  the  analysis  of  the  coupling  of  a  single  collision  com¬ 
plex  to  radiation  fields  in  the  Hilbert  space  of  the  scattering  states..  The  details  of  vari¬ 
ous  derivations  and  further  discussions  of  the  implications  of  the  scattering  state  repre¬ 
sentation  are  given  in  Ref.  15. 

The  essence  of  the  new  method  is  to  pair  individual  atoms  such  that  one  can  formally 
solve  for  the  binary  collisions  by  describing  the  quantum  state  of  a  pair  by  the  scattering 
state  corresponding  to  a  specified  asymptotic  state  of  the  pair  in  the  limit  of  vanishing 
collision  potential.  If  HA  represents  the  asymptotic  Hamiltonian  of  the  two  atoms  that 
are  paired,  and  V  is  the  localized  collision  potential  [Ref.  16],  then  the  scattering  state 
corresponding  to  an  asymptotic  state  \rp)  is  given  by  [Refs.  16-18] 

|0)+  =  n(+)iv>)  (1) 

n(±)=l  +  {E -HA  -V  ±ie)~1V  (2) 

where  E  is  the  energy  of  the  asymptotic  state  j xji)  and  is  treated  as  a  c- number  param¬ 
eter.  The  transformation  Equation  1  can  be  viewed  as  a  canonical  transformation  when 
there  are  no  bound  states  of  the  two  atoms  induced  by  V,  or  even  if  there  are  such  bound 
states  for  an  attractive  V,  provided  they  can  be  neglected  with  little  effect  on  the  dy¬ 
namics  of  pairs.  This  follows  from  a  fundamental  theorem  [Ref.  16]  of  the  Quantum 
Scattering  Theory: 
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Theorem:  To  every  solution  x />t  of  the  Schrodinger  equation 


ihj^  =  (Ha  +  V)xl>t  (3) 

which  is  orthogonal  to  all  bound  states,  there  corresponds  a  unique  solution  <$>t  of  the 
particle  Schrodinger  equation 

(4) 

(5) 

Conversely,  to  every  solution  of  Equation  4  there  corresponds  exactly  one  solution  of 
Equation  3  such  that  Equation  5  holds. 

The  theorem  shows  that  the  asymptotic  pair  states  can  be  put  into  one-to-one  corre¬ 
spondence  with  the  scattering  states  in  the  absence  of  bound  states: 

l*.(l)>®!*(2)>«— |*(1,2)>+  (6) 

Here  the  asymptotic  states  of  the  pair  are  written  as  direct  products  of  atomic  states, 
since  the  two  atoms  evolve  freely  when  the  collision  potential  vanishes.  In  fact,  fI(  +  )  sat¬ 
isfies  [Ref.  17j 

(n(+))tn(+)  =  i  (7) 

n'+'(nl+l)'  =  i-£>,W,|  (8) 

$ 

where  |t^)’s  are  the  bound  states  of  HA  -f  V .  Thus,  fl(  +  )  is  actually  a  unitary  operator 
and  Equation  I  becomes  a  unitary  transformation,  if  the  bound  states  are  ignored. 


such  that  the  norm 


<j>t 


lim  IJ0t  -  (j>t  |!  =  0 
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If  the  bound  states  are  to  be  neglected,  then  there  must  be  no  chemical  phase  transitions 
in  the  gas.  In  other  words,  the  temperature  and  the  pressure  of  the  gas,  as  well  as  the 
type  of  atoms  (or  molecules),  must  be  such  that  no  chemical  transformations  take  place 
locally  or  otherwise.  The  gas  may,  of  course,  be  a  mixture  of  different  chemical  species. 
But  each  species  in  the  mixture  must  preserve  its  chemical  identity,  and  there  must  be 
no  phase  transformations  from  one  species  to  another. 

In  the  asymptotic  limit,  the  operators  describing  various  properties  of  a  collision  com¬ 
plex  can  be  obtained  from  the  atomic  operators.  To  go  to  the  scattering  state  represen¬ 
tation,  the  asymptotic  collision  complex  operators  are  transformed  by  n(  +  ),s  and  be¬ 
come  energy  and  potential-dependent.  For  example,  let  A  be  an  atomic  operator  such 
that  the  corre5ponding  operator  of  the  complex  is  additive,  as  given  by  A(l)  +  A{2).  In 
the  scattering  state  representation,  A  is  transformed  into  + 

A-+  A(  +  )  =[1  +  V{E'  ~Ha  -V  -*£)-1]A[1  +  (£-Fa  —  V  +  te)~ 1 V]  (9) 

If  A  is  some  moment  operator  (electric  dipole,  magnetic  dipole,  etc.),  then  A(+)  couples 
to  radiation  fields  in  the  scattering  state  representation.  In  contrast,  in  the  standard 
theories,  it  is  A  that  is  coupled  to  radiation  fields.  The  form  of  Equation  9  shows  that  the 
ordinary  symmetry  arguments  and  selection  rules  of  the  standard  theories  based  on  A 
do  noi.  hold  in  the  present  approach.  The  scattering  state  representation  mixes  the  sym¬ 
metry  properties  of  the  pair  of  asymptotic  atomic  states  with  the  collision  potential.  In 
general,  y4(  +  )  has  mixed  symmetry  with  respect  to  rotations,  inversions,  etc.,  even  if  A 
has  a  definite  symmetry. 

In  the  Kroll- Watson  theory  of  inelastic  atom-atom  scattering  in  a  laser  beam  [Ref.  19], 
the  coupling  of  the  collision  complex  to  the  radiation  field  is  analyzed  by  expanding  the 
state  of  the  complex  in  terms  of  the  actual  molecular  states,  with  time-dependent  ampli¬ 
tudes.  These  molecular  states  are  the  bound  states  of  the  two  colliding  atoms  with  po¬ 
tential  V.  In  the  present  theory  these  bound  states  are  completely  excluded.  The  quasi- 
molecular  operators  are  obtained  from  the  free  atomic  operators  by  transformations  of 
the  type  in  Equation  9. 
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Consider  the  Hamiltonian  for  two  indistinguishable  atoms  colliding  in  the  presence  of  ra¬ 
diation  fields.  The  assumption  of  indistinguishability  is  not  necessary  for  the  formalism; 
but  it  simplifies  most  of  the  results.  This  assumption  cam  be  readily  removed  and  the  re¬ 
sults  adapted  to  the  case  of  two  nonidentical  atoms.  When  the  center  of  mass  (c.o.m.) 
motions  of  the  atoms  are  separated  from  their  internal  motions,  the  asymptotic  state  of 
the  pair  can  be  written  as 


—  {XiX2\  AiA2\ciiPiCl2p2o) 

+  oxl>„,  (A,  )11„ ,  (A,  )«‘ < ■*'  * f'  * ' 1 1 


(10) 


where  X1<2  are  the  atomic  c.o.m.  positions,  pii2  are  the  atomic  c.o.m.  momenta  (di¬ 
vided  by  h),  A1i2  are  the  atomic  internal  coordinates,  and  ai,2  are  the  atomic  internal 
quantum  numbers.  Vol  is  the  volume  to  which  the  pair  is  confined.  Since  the  atoms  are 
indistinguishable,  the  Hamiltonian  is  symmetric  with  respect  to  the  interchange  of  the 
two  atoms.  It  follows  that  the  asymptotic  eigenstates  of  the  complex  can  be  symmetric 
or  antisymmetric  with  respect  to  the  interchange  of  the  two  atoms.  To  allow  for  both 
possibilities,  the  variable  a  is  used.  It  can  assume  only  two  values,  +1  and  —1.  The 
c.o.m.  quantum  numbers  pi  and  p2  are  used  interchangeably  with  the  total  momentum 
of  the  pair,  p  =  pi  +  p2 ,  and  the  relative  momentum  of  the  atoms,  q  =  pi  —  p2 .  The 
designations  for  the  pair  states  |ai pi  a2p2a)  and  lo^  ol2o\ pq)  are  used  interchangeably. 
The  collision  potential  V  depends  on  the  atomic  internal  oordinates  Ax  and  A2 ,  and  on 
the  interatomic  separation  Xt  —  X2 .  V  is  symmetric  with  respect  to  the  interchange  of 
the  atoms.  Its  Fourier  matrix  elements  are  given  by 


V. 


^  K,  J 


A  i  A  : 


dRe' 


k-R 


X  rai  Ml  w;,  {A2)V{A,A2]R)^  (a3) 


(11) 
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The  coupling  of  the  complex  to  radiation  fields  is  described  by  multipolar  couplings. 
Specifically,  the  dipole  moment  density  is  written  as 


=  I>(r-  *y)<5  (12) 

3=1 

where  d,  is  the  electric  dipole  moment  operator  of  the  jth  atom.  Although  this  report 
focuses  on  radiative  couplings  that  are  given  exclusively  in  terms  of  the  electric  dipole 
moments,  the  results  can  readily  be  generalized  to  magnetic  dipoles,  quadrupole  mo¬ 
ments,  etc.,  by  adding  further  terms  to  Equation  12  and  by  treating  D  as  a  generalized 
moment  density.  The  Fourier  transform  of  D  is  given  by 

I  dre'*  *  D(f)  =  =  D{k)  (13) 

J  3=1 

The  asymptotic  energy  eigenvalues  of  the  complex  axe  designated  by  J5ai?l0l p, : 

-fiMaiPi«2P»  =  Ea^ia7tFi\a1pla2p2a)  (14) 


(15) 


where  Eai  and  Eai  are  the  atomic  internal  energies.  The  energy  of  the  complex  has  the 
same  form  (Equation  15)  regardless  of  the  sign  of  a.  The  scattering  state  corresponding 
to  IctiP! o-2p2a)  is  designated  with  a  plus  sign: 


\diPiQ.2p2a)Jr  =  1 +  (£«,?,«, ?, 


A 


aiPia2p2a) 


(16) 


Using  the  completeness  of  the  scattering  states  (in  the  absence  of  the  bound  states),  the 
total  Hamiltonian  is  written 


H  =  Ha  +V  +  H,+Hint 


(17) 
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where  H ,  is  the  Hamiltonian  for  the  free  radiation  fields,  as 


H=  ^2  |<*iPi«aft*}+  +  («iPi«afi*l 

<*  101“  lPl* 

+  >  ku.za?  _a.  r 

fa  ^  AJfc  ^  ^ 

\k 

+  |*  5Z  (PiP^P/iPa') 

x  |aipiot2ftff)+  +  (aiPi'a^Pa^l  +  herm.  conj.j  (18) 

where  axi  and  a*  -  are  the  annihilation  and  creation  operators  for  the  mode  propagating 
with  momentum  k,  the  polarization  index  A  •  exi  is  the  polarization  vector,  and  u>x-  is 
the  frequency  of  the  mode.  The  radiative  coupling  coefficients  are  defined  by 

A  f  ,  „  (2irhux-k\> 

,  aj  ja'j  a  j  vP1  i  Pi  P2  /  (  y  J 

x  +  {a1p1a2p2a\exS  •  D[k)\a'lpxa'^p:;a)+  (19) 

The  matrix  elements  of  D  between  two  asymptotic  states  of  the  complex  are  given  by 


( <* i p ! p2 a | P ( fc) | a; p2' a)  =  $(Pi +  k)6[p2,p±)da^ 


+  ,  a (pi  ,  Pi'  )  6  (p2  ,  PA'  +  k)  da  ,  a  ; 


+  o 


5a,a'/(Pi,P2  +  k)6{p2ip')daia't 


+  6aia'J(Pi,  P2)S(p2 ,  Pi  +i)4,.[ 
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=  <5(P,P'  +  £)  j^.a; 6(9,4  '  +  fc)<Taia.  +6aia't6(q,q' '  -  £)<£,, 


+  a 


(20) 


The  subscripts  on  d,  are  omitted  here,  because  the  atoms  are  indistinguishable,  and  the 
atomic  matrix  elements  are  designated  as  (ajdja')  =  daa>.  The  matrix  element  in  Equa¬ 
tion  20  represents  a  transition  in  which  a  photon  is  absorbed  and  the  total  momentum  of 

the  complex  is  increased  by  k.  The  first  two  terms  are  a  pair  of  delta  functions  located 
— 

at  q  =  q‘±k.  They  describe  photon  absorption  by  one  or  the  other  colliding  atom. 

The  other  two  delta  functions  at  q  =  —q'  ±  k  arise  from  the  absorption  of  the  pho¬ 
ton  accompanied  by  the  exchange  of  the  atoms.  These  singularities  may  be  called  the 
exchange  singularities.  If  the  atoms  were  not  paired  and  just  one  atom  was  being  con¬ 
sidered,  there  would  be  just  one  delta  function  singularity,  for  example  at  px  =  pi  '  +  k. 
Pairing  of  the  atoms  doubles  the  number  of  singularities.  If  the  paired  atoms  are  iden¬ 
tical,  the  number  of  singularities  is  doubled  again.  If  the  photon  momentum  is  small 
compared  to  the  relative  momenta  of  the  atoms,  \k\  <C  |g  |  cz  \q  '|,  then  the  two  singu- 
laxities  at  q  =  q  '  ±  k,  as  well  as  the  exchange  singularities,  merge  together.  Note  that  |A;| 
may  be  completely  negligible  compared  to  the  total  momentum  of  the  complex,  yet  be 
significant  compared  to  the  relative  momentum  of  the  atoms.  Merging  of  the  singulari¬ 
ties  occurs  only  if  |fc|  is  also  negligible  compared  to  the  relative  momentum.  Finally,  for 
a  given  set  of  quantum  numbers  <xx  a2  ot\  oc'2 ,  the  actual  number  of  singularities  appear¬ 
ing  in  the  asymptotic  matrix  element  of  D  depends  on  how  many  of  the  atomic  matrix 
elements  are  nonzero. 

Collisions  modify  the  preceding  physical  picture  significantly.  The  matrix  elements  of  D 
in  the  scattering  state  representation  are  given  by 
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+  {alp1a2p2c\D{k)\a[p{a!2p2a)+  =  {a^  a2P2<r\D\a[p1  'a2p2  'a) 


+  {aiPiCt2p2o\V (Eaifxatfi  -Ha  -  V  -  it)  V\a\pt  '  a!2p2  ' a ) 

+  {^P^iP^<y\v[Ea‘xfx  ,a-j<  -Ha  -  V  +  it)  V\ a[p'o^p2a) 

+  (<*iPi<*2P2<r\V  (Eaxifxatfi,  —  Ha  —  V  —  it)  V 

x  [E<K<K  ~  Ba  -  V  +  it)~lV\a\pla'2p'a).  (21) 

The  collisional  parts  of  Equation  21  cam  be  evaluated  by  inserting  the  complete  set  of 
the  asymptotic  states  at  appropriate  places  and  using  Equation  20.  To  lowest  order  in 
V,  the  collisional  part  of  Equation  36  is  given  by 

{diPia2p2a\6V[k)\a\ p' a2 p2 a)  ce  j 6{p,p‘  +  k)^ 
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A  comparison  of  Equation  22  with  Equation  20  shows  that:  (l)  The  number  of  singular¬ 
ities  is  increased:  in  other  words,  new  coupling  channels  are  opened  up  by  the  scattering 
process.  The  new  singularities  are  simple  poles,  arising  from  the  denominators  in  Equa¬ 
tion  22.  They  form  surfaces  in  the  space  c  the  ’^tive  momentum.  If  the  Fourier  trans¬ 
form  of  the  collision  potential  has  singularities,  then  these  also  contribute  to  the  cou¬ 
pling  coefficients.  (2)  The  original  as  well  as  the  new  singularities  are  broadened.  This 
broadening  is  due  to  the  appearance  of  simple  pole  singularities.  The  Fourier  transform 
of  the  collision  potential  also  contributes  to  this  broadening. 

Of  course,  collisionally-opened-up  radiative  transitions  have  been  observed  for  a  long 
time  [Ref.  20]  and  are  well-known  from  collision-induced  absorption  [Ref.  12],  field  in¬ 
duced  avoided  crossing  [Ref.  13],  pressure-induced  four  wave  mixing  [Ref.  10],  etc.  The 
difference  between  the  present  theory  and  others  is  that  in  the  other  formalisms  V  en¬ 
ters  into  the  physical  picture  indirectly,  either  through  linewidths  and  lineshifts  [Ref.  10] 
or  through  the  calculation  of  a  particular  molecular  state  wave  function.  In  the  present 
formalism,  V  is  directly  incorporated  into  the  coupling  coefficients  of  atoms  to  radiation 
fields.  These  coupling  coefficients  take  into  account  the  collisions  between  the  atoms  of 
the  complex  to  all  orders  in  V . 

A  further  discussion  of  the  coupling  coefficients  between  the  complex  and  radiation  fields 
is  given  in  Ref.  15. 

In  order  to  obtain  an  analytically  solvable  Hamiltonian,  some  simplifying  assumptions 
are  required.  Assume  that  the  atoms  of  the  complex  are  two-level-atoms  and  that  the 
complex  is  coupled  to  just  one  radiation  mode.  Furthermore,  assume  that  the  momen¬ 
tum  exchanged  during  collisions,  as  well  as  the  photon  momentum,  is  negligible.  The 
latter  assumption  means  that  the  atoms  of  the  complex  preserve  their  individual  c.o.m. 
momenta  at  all  times,  which  is  essentially  equivalent  to  the  straightline  approximation  of 
many  of  the  semiclassical  theories  of  atomic  collisions.  Finally,  assume  that  a  =  +1  and 
omit  references  to  a  in  the  following  discussion. 

From  Equation  10,  there  are  four  asymptotic  states,  which  can  be  written  as 
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10)  =  ^[|0;pi)A  ®  |0;p2)B  +  |0;pi)B  ®  |0;p2)A]  (23) 

11)  =  [l°; Pi )a  ®  |l;p2)s  +  (0; pi ) a  ®  |l;p2)>t]  (24) 

1 2)  =  ^[)1;Pi)a  ®10;p2)a  +|l;Pi)fl  ®  |0;p2)A]  (25) 

13)  =  4|[|1;p*)a  0  11;p*2)b  +  |i;p*i)b  ®  |i;p2)A]  (26) 


Here,  the  state  vectors  |o:;pI)A,a  (a  =  0,1;*  =  1,2)  refer  to  the  individual  atomic  states. 
a  designates  the  internal  atomic  state,  p<  the  c.o.m.  momentum.  A,  B  label  the  atoms. 
One  should  interpret  the  meaning  of,  for  example,  lOjpj)*  as 

|0,Pi)A  —  (Aa ; XA  (0; px )  =  '^^(4.)  (27) 

VVol 

where  Aa  b  refers  to  the  sets  of  the  internal  variables  of  the  atoms  A  and  B,  and  XA  B 
to  their  c.o.m.  coordinates.  ip0  is  the  wave  function  for  the  internal  ground  state.  Un¬ 
der  the  present  approximations,  pi  and  p2  axe  constant  parameters  without  dynam¬ 
ical  significance.  It  is  therefore  entirely  sufficient  to  label  these  asymptotic  states  by 
—  0, 1,2,3).  They  are  the  eigenstates  of  the  asymptotic  Hamiltonian  HA  : 

»<W=P,W  (28) 

where 


El  =  2  E0  + 


h2p[  ,  ft3?! 


2m 


+ 


2m 


h2pl  h  2p2 

E\  =  El  =  E0  +  Ei  +  +  —  — 


2m 


2m 


(29) 


(30) 
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h2Fi  ,  h2Fi 


El  =  2  El  +  -^  + 


2m  2m 


(31) 


The  corresponding  scattering  states  are  given  by 


!**>♦  =  [i  +  (b;  -  ha  -  v  +  «)-  ‘  v]  |„> 


(32) 


which  are  the  eigenstates  of  HA  -f  V : 


+v)\h)+  =  £»+ 


The  simplified  Hamiltonian  obtained  from  Equation  18  can  therefore  be  written  as 


(33) 


H  =  ^  E"  \n)+  +  (vt\  +  hua+  a  + 


|f)++  (o| 


+  920 12)++  (0|  +  ^3i  |3)+  +  (l|  +  ?32 13)+  +  (2|)d  +  herm.  conj. 


(34) 


where 


9„»'  =  (^)'  +  M«-PVK  (35) 

In  writing  Equation  34,  the  rotating  wave  approximation  is  made  and  the  coupling  be¬ 
tween  |1)+  and  |2)+  is  omitted,  since  they  have  the  same  energy. 

The  eigenvalues  and  eigenstates  of  the  Hamiltonian  in  Equation  34  can  be  obtained  by 
writing  an  arbitrary  state  \rp)  as  a  linear  superposition  of  the  states  of  the  complex  and 
the  radiation  field,  in  the  form 


NO  =  lo  |0)+  ®  jn  +  2)  +  |l)+  ®  |n  +  1)  +l2\2)+  ®  |n  +  1)  + /3 13)+  ®  |n) 


=  /0|«0)+  +  li  |nl)+  +  /3|n2)+  +  /3  |n3)+ 


(36) 
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where  |n)  designates  the  radiation  state  with  n  quanta  (n  >  0).  From 


H\tJ>n)  =  E\t(>n)  (37) 


one  finds  that 


Eq  +  (n  +  2)tiu 

—ihg\  0y/n  +  2 

—th&  0y/n  +  2 

0 

Iq 

lo 

ihgu>Vn  +  2 

El  +  (n  +  l)^w 

0 

-ihg^y/n  +  1 

h 

=  E 

h 

ihg30y/n  +  2 

0 

E^  +  (n  ■+■  l)  tiui 

-ihgl2y/n  +  1 

l* 

h 

0 

ihg3ly/n+  1 

ihgaiVn  + 1 

El  +  nhuj 

h 

h 

(38) 

Thus,  the  eigenvalue  problem  for  H  is  reduced  to  the  diagonalization  of  a  4  x  4  hermitian 
matrix.  Denote  this  matrix  by  Hn .  From 


det(tf„  -E)  =  0 


(39) 


and  the  definitions 


E  =  h£  +  El  +  (n  +  l)ftu; 


(40) 


h&  =  El  -  Ec0  -  huj 


(41) 


one  finds  the  eigenvalue  equation 
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o  =  ?  -  e 

-£A 


A3  +  (n  +  1)(|^3X  T  +  1^32  I3)  +  (n  +  2)  (1^10  |3  +  1^30  r) 
(n  +  l)(|03i|3  +I&32I2)  -  (n  +  2) {\gio  |3  +  |02o  |a ) 


+  (n  +  l)(n  4-  2} 


1010  |3  |033  |3  +  1^20  |3  |ff31  |3 


—  010  020031  033  —  010  020031033 


(42) 


The  indistinguishability  of  the  atoms  of  the  complex  and  the  straightline  trajectory  as¬ 
sumption  introduce  certain  symmetry  relations  among  which  simplify  Equation  42. 
These  are  [see  Ref.  15]: 


Similarly, 


001  —  002 


(43) 


013  ~  023  (44) 

It  follows  from  Equation  43  and  Equation  44  that  the  last  term  of  the  eigenvalue  equa¬ 
tion  vanishes.  One  £  factors  out.  Designate  the  roots  of  Equation  42  by  f„A,  where 
A  =  0, 1,2,3.  The  £  -  factor  gives  a  root  that  is  zero.  Let  this  root  be  f„3.  The  remain¬ 
ing  three  roots  are  the  solutions  of  the  cubic  equation 


0  —  £3  —  [2(n  +  2)(<7io  |3  +  2(n  +  1)(<73i  |3  +  A3]  £ 

+  2A[(n  +  2)(<710|3  -  (n  +  l)(03i  |2]  (45) 

Only  the  real  solutions  of  Equation  45  are  acceptable  for  energy  eigenvalues.  The  num¬ 
ber  of  real  solutions  of  Equation  45  depends  on  the  sign  of  the  quantity 
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3  2 

tn  =  U  —  V 


where 


tt„  =  -[A3  +  2(n  +  2)|g10|2  +2(n  +  1)|<73i|2] 


vn  =  A[(n  +  2)|^10|3  -  (n  +  1)|$31|3]  (48) 

Since  Hn  is  hermitian,  wn  >  0,  the  cubic  equation  is  irreducible  and  there  axe  three  real 
roots  [Ref.  20]: 


Zno  =  2^/tCcos 
Cm  =  — 2\/tt^  cos  (^-  g—  ^ 
Cr»2  =  -2,^ COS  ^  ^ 


where 


=  fan' 


■(=?) 


Let  £„au  represent  the  components  of  the  eigenvector  corresponding  to  the  eigenvalue 
Cn a  •  If  the  eigenvalues  are  not  degenerate,  then 
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(53) 


tfiXO  —  T7  [£n\  (fnA  —  A)  —  2  (n  +  l)|<?31  |2  ] 

*'n\ 

ln\i  =  lnx2  —  Tp — Vn  +  2(£„a  —  A)ffio  (54) 

•‘’nA 

2  - 

in  A  3  = 77"  —  V  (j1  +  l)(n  +  2)g10<73i  (55) 

■i’nA 


where 


JV„a  =  {[&*(£»*  -A)-2(n  +  l)|ff3ir]3+2(n  +  2)(ef,A  -A)3|g10|3 

N  1/2 

+  4(n  +  1) (n  +  2) |y10 13  \g3i  |3  j 

If  the  eigenvalues  are  degenerate,  an  orthonormal  set  of  eigenvectors  may  be  determined 
directly  from  the  matrix  Equation  38.  Because  H„  is  hermitian,  the  eigenvectors  are 
complete  as  well  as  orthonormal: 


a*  (57) 


J2lnxJnXu'=^>  (58) 


The  eigenstates  of  H  are  given  by 


|0„a)  =  ]T]/„aJ  nfi)+ 


(59) 


H\^)nx)  ~  Enx\lpnx)  ~  [hC„A  +  E{  +  (n  +  l)ftu/]  |0„A  )  (60) 
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An  arbitrary  state  of  the  collision  complex  plus  the  radiation  field  can  be  written  as 

w  =  (61) 

n  =  0  A  =  0 

If  the  state  of  the  total  system  at  time  t  =  0  is  given  by  a  state  like  Equation  61,  then 
the  state  at  time  t  is  given  by 


WO)  W  =  (62) 

nX 

The  corresponding  density  matrix  operator  is  given  by 


p{t)  =  |V>(*))W)I  =  J2 


c<  +  )  (+)< 
CnA  Cn'A' 


n'A’ 


(63) 


n  A  n '  A ' 


If  at  time  t  =  0  the  system  can  be  described  only  by  a  mixture  of  states,  for  example  by 
the  density  matrix  operator  p(0),  then 


0(0)  =  £ 


'X' 


nAn'V 


(64) 


o(t)  =  £  !*•»)[«  <( 


^a!«*a*(°)](^'a* 


nAn'A' 


nAn'A' 


(65) 


Clearly,  the  specification  of  either  the  initial  state  amplitudes  c^1  in  Equation  62  and 
Equation  63  or  the  initial  density  matrix  elements  P^!n.A.(0)  in  Equation  64  is  equiva¬ 
lent  to  a  complete  solution  for  the  dynamics  of  the  system. 
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The  initial/boundary  conditions  are  specified  in  the  asymptotic  limit  and  for  individ¬ 
ual  atoms.  It  is  assumed  that,  at  t  =  0,  the  density  matrix  of  the  total  system  can  be 
factored  into  an  outer  product  of  the  density  matrix  for  the  complex  and  the  density 
matrix  for  the  field.  These  assumptions  correspond  to  a  physical  picture  in  which  the 
radiation  field  and  the  collision  potential  are  turned  on  simultaneously  at  t  —  0.  From 
Equations  23-26 

I P>  =  Calix;a,i,\al'Pii  >*  ®\<*2'>Pi>  >B  (6« 

a  1*1  <X3«3 

where  a  —  0,  1  designates  the  internal  quantum  states  of  the  atoms  and  *  =  1,  2  desig¬ 
nates  the  two  momenta  for  the  atoms.  All  of  the  nonvanishing  C^\.aV  ’s  are  equal  to 


^03;01 


r3 

^11:12 


_  _  /->1  _  (~<2  _ 

~  '-'01:12  —  v"/12;01  —  ^11:02  —  °02;ll 

-  c3  —  — 

"0l2;U-v/2 


From  Equation  59, 


Since  both  the  collision  potential  and  the  radiation  field  are  turned  on  at  t  =  0, 


^r»A;n'A’ (0)  ~  AM  ^n'X  V'  P„  +  m  *  ,n'  +  m  ,  (®) 
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P^'(°}  -  S  Pt, (°)p«, i, 


(70) 


“  1  1 1 “1 1  J 


Here  p7  (0)  is  the  density  matrix  operator  for  the  field  at  t  =  0,  p°  (0)  is  the  density  ma¬ 
trix  operator  for  a  single  atom,  and  is  an  integer  defined  by 


v  =  { 


2  for  n  =  0 

m„  =  <  1  for  n  =  1,2 

0  for  p  —  3 


(71) 


If  the  radiation  mode  is  in  a  coherent  state  and  the  average  number  of  quanta  is  n,  then 


pI„.(o)  = 


«-*(a) 


(72) 


Assume  that  pa  (0)  is  diagonal  and  the  occupation  probabilities  are  given  by 


p:,;,.,.(o) = 6,„'S,.’h  (p.) 


(73) 


where  fa  (p*)  is  the  conditional  probability  that  the  internal  state  of  the  atom  is  a,  if  the 
atom  is  in  the  c.o.m.  momentum  state  pf.  For  these  initial/boundary  conditions 


MM 


a-f  m  «  +  n#  +  * 


y/{n  +  m„)\(ri  +  mM-)! 


X  InXuln'X'n'FnH'iPlyfa) 


(74) 
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where 


=  /o(p,)/o(p1)n‘:"(E;)ni;'(E;,) 

+  /o(p-.)/.(ftK:”(£;)n';>(£;.) 

+  A  (ft )  A  (ft )  n£  > 1 ’  (e;  )  n'  *  ’  (e;  . ) 

+  (75) 

Taking  the  trace  over  the  radiation  field  in  Equation  74  and  75  and  transforming  back  to 
the  Hilbert  space  of  the  asymptotic  states, 


(76) 


where 


P*t'(0  =  «"*  53  *(n-f  nv,n'  +  nv) 


x  e 


nn'U'p'V" 


(n) 


y/{n  +  mM-.)!(n'  +  mM«.)! 


x  *r»  A*  C  A  ‘»'V#*"'C’A' 


(77) 


The  elements  of  the  atomic  density  matrix  operator  are  given  by 


Pai-.a'i'i*)  —  53  W  (^aV;a"i") 


(78) 


The  physically  interesting  case  is  the  one  in  which  n  »  1  and  the  probability  distribu¬ 
tion  of  n’s  in  Equation  72  is  sharply  peaked  around  h.  If  the  presence  of  a  moderately 
intense  coherent  field  is  assumed,  and  also 


(79) 
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(80) 


Vurt.  =  -v;„ 


00 


to  first  order  in  V  and  in  the  limit  €  — ►  0: 


•““®*|-<*(1  +  5T^)+0(V*) 


(81) 


(82) 


It  follows  that,  to  first  order  in  V : 


93i  -  9io  =  9 


(83) 


Under  these  assumptions,  the  eigenvalue  equation  becomes 


e4-(A2  +  4%i2)e2  =  0 


(84) 


The  eigenvalues  are 


£o,i  =  ±y/  A2  +  4n|<j|2  = 


(85) 


$3,3  ~  0 


(86) 


i 

j 

The  eigenvectors  corresponding  to  are,  respectively,  ! 
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<'-A 

2f* 

2«' 

frgy/2n 

Ji  = 

j7gy/2fi 

frg\/2n, 

jrgV 2i 

_  <'+* 

t'+± 

3  (' 

2 1' 

The  eigenvectors  corresponding  to  0  are 


-  pg'Vzn 

0 

A 

i 

h  * 

y/2 

A 

2 

_  l 
\/  2 

-  jrff* 

0 

.  « 

Under  the  present  approximations,  Equation  76  becomes 


(87) 


(88) 


AA' 


(89) 


which  can  be  rewritten  as 


/C'(0  =e 


_  m  .  )w« 


-  1 


+  R'^.e 


2»{  t 


-2 


+  ^(.'e 


(90) 


The  explicit  expression  for  ( t )  is  given  in  Ref.  15,  where  Equation  (4.48)  and  Equa¬ 

tion  (4.49b)  show  that  the  quantity  is  the  generalized  Rabi  frequency  for  an  atom  in 
the  coherent  fundamental  mode  [Ref.  22j.  If  8  is  the  angle  between  the  polarization  vec¬ 
tor  and  the  bare  atomic  dipole  moment 
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€•  dl0  =  |dlo|cOS0 


(91) 


then  £'  can  be  written  as 


C  =  |A|  v^l  +  ?2  cos2  6 


(92) 


where 


fa  =  f^ls(*=0)|,  =  2s-  (93) 

is  the  square  of  the  ratio  of  the  standard  Rabi  frequency  to  detuning,  whose  magnitude 
is  modified  by  the  collision  process.  It  is  seen  from  Equations  81-83  that  the  effect  of 
the  collision,  to  first  order  in  V,  is  to  modify  the  magnitude  and  the  phase  of  the  atomic 
dipole  moment. 

The  appearance  of  just  one  angle  variable  in  Equations  91-93  may  be  surprising  at  first 
glance,  since  there  are  two  atoms  in  the  complex.  However,  it  was  assumed  that  the  two 
atoms  of  the  complex  were  indistinguishable.  If  the  internal  d10  were  to  point  in  differ¬ 
ent  directions  for  each  atom,  then  the  atoms  would  be  distinguishable.  The  symmetric 
and  antisymmetric  combinations  of  asymptotic  atomic  states  can  still  be  used  for  distin¬ 
guishable  atoms.  For  symmetric  combinations,  the  replacement  is 

?  •  <T10  — ♦  -e  •  (dl[A)  +  <T1(oB))  (94) 

in  Equations  81-83  for  the  distinguishable  atoms  A  and  B.  Formulas  for  distinguishable 
atoms  are  more  complicated  than  those  for  indistinguishable  atoms.  The  matter  is  not 
pursued  in  this  report. 
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The  initial  conditions  are  incorporated  into  the  density  matrix  solutions  through  If 
the  matrix  elements  of  nl  +  )  are  expanded  in  powers  of  V,  then 


Fftfi  —  ft o  fo  (Pi )  /o  (Pa )  +  ft  i  fo  (Pi )  fi  (P2 ) 

+  fta  fo  (ft  )/i  (ft )  +  ft3  h  (ft  )/i  (ft )  +  0(F2 )  (95) 

In  other  words,  up  to  and  including  the  first  order  in  V,  is  given  by  the  conditional 
probabilities  fa  (ft).  The  latter  is  readily  obtained  when  the  atoms  of  the  gas  obey  ei¬ 
ther  the  Boltzmann  statistics  or  the  Bose-Einstein  statistics  (since  it  was  assumed  that 
a  —  +1  with  respect  to  the  interchange  of  the  atoms  of  the  complex).  For  the  Boltz¬ 
mann  statistics,  the  probability  is  that  an  atom  has  the  c.o.m.  momentum  ft  and  that  the 
internal  quantum  state  a  is  proportional  to  exp  (3(Ea  +  ^£J-)| ,  where  (3  =  (kB  T)~ 1 . 
Consequently,  the  conditional  probability  that  the  atom  with  ft  is  in  the  state  a  is  given 
by 


/.(*)« 


(96) 


The  /„ ’s  are  thus  independent  of  the  c.o.m.  momenta.  Hence, 


Foo=fl ,  fu=f»=/oA,  Fzz  =  fl 


(97) 


The  fa ’s  are  independent  of  ft  for  the  Bose-Einstein  statistics  as  well: 


fa 


1 


Ef  «**«-!  N 

•'  \ «**«*'  -  1  / 


(98) 


The  spectrum  of  the  diagonal  matrix  elements  of  pe(t )  consists  of  five  peaks  at  0,  ± 

±2  The  peaks  at  ±  £'  and  ±  2£'  disappear  if  f  C  1,  that  is,  if  the  frequency  detun¬ 
ing  from  the  resonance  is  large  compared  to  the  Rabi  frequency: 


f<l  => 


Pit  (0  —  fl , 

p£(0— -MWo), 


pSJW— ‘(Wo)2  (99) 

In  other  words,  /?'+  (t)  is  not  changed  from  p^+  (0),  because  the  coupling  to  the  funda¬ 
mental  mode  is  negligible.  On  the  other  hand,  near  or  at  resonance,  p„£  and  have 
five  peaks,  but  p\X  and  pc22  do  not  have  the  components  corresponding  to  ±£': 

i=> 


Poo  (0  — *  \  +  £  (/o  ~  |)  +  (/o  -  cos(flR  |  cosfljt) 
+  ^  (/o  -  ^  cos(2n,  |  cos  0\ t) 

pstw./£w  — 

-  ^(/°  ~  0  cos(n*|cos0|f) 
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Pit  (Oi  — *  J  +  \  (/o  -  -  (jo  -  0  cos(nfl  |cos0|t) 

+  \  (/o  -  ^  cos(2 n*  I  cos  0|t)  (100) 

The  limits  S'  — ►  oo  and  f  — *  0  will  be  referred  to  as  the  strong  and  the  weak  coupling 
regimes,  respectively. 

The  coherent  radiation  field  creates  off-diagonal  correlations  in  the  density  matrix  of  the 
complex.  The  spectral  lines  of  p%+ ,  peQ+ ,  p'+ ,  and  p\+  are  at  w,  ±  and  w  ±  2£\  The 
lines  of  p«+  are  at  0,±£\  and  ±2f .  />;£  has  spectral  lines  at  2w,  2w  ±  and  2w  ±  2£\ 
These  correlations  arise  from  the  coupling  of  the  complex  to  a  coherent  field.  The  role  of 
the  collision  potential  is  to  modify  the  coupling  coefficients  to  the  radiation.  In  the  weak 
coupling  limit,  all  of  these  off-diagonal  matrix  elements  vanish:  pc0  +  and  p\+  as  f2 ;  the 
others  as  £.  In  the  strong  coupling  limit: 

£—►00  => 


(/o-i) 

+  (/,  —  -)  sin{2fifl  |  cos  S\t) 


sin(flR|cos0|t) 
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Pit  —  \*2iut  (/o  -  [l  -  cos(2nfl  \cose\t)]  (101) 

where  g  =  \g\  exp(t0„).  The  first  four  elements  in  Equation  101  oscillate  at  w  ±  ClR  |  cos  6 1 
and  u)  ±  2ClR  |  cosfl|,  as  could  be  expected  from  single  photon  transitions  and  from  the 
fact  that  both  of  the  atoms  are  executing  Rabi  oscillations  with  frequency  ClR  |  cos  0|.  It 
is  interesting,  however,  that  for  sufficiently  large  f ,  which  corresponds  to  polariza¬ 
tion  at  the  second  harmonic,  oscillates  at  2ui  and  2ui  ±  2flfl  |  cos 0|,  the  components  at 
2 u)  ±  fin  |  cos  9\  are  suppressed,  and  that  p\%  becomes  independent  of  the  radiation  field. 

Because  of  the  form  of  fi(+)  pa(t)  can  be  decomposed  into  tliree  parts.  From 

Equation  2, 


1 


+ 


\v 


Eh'"  —  Ep‘1  —  V  —  if 


\p")  +  1 


Eu.-Eu-V  +  ie 


rV\p.') 


+  <H 


Ep‘  —  E^  —  V  +  if 


-v\/iW \v 


E^i‘1  —  E,t"  —  V  —  te 


r>") 


(102) 


Here  the  terms  are  grouped  according  to  the  number  of  kronecker  deltas.  The  term  with 
two  kronecker  deltas  is  zeroth  order  in  V.  The  terms  with  one  kronecker  delta  are  at 
least  of  the  order  of  V.  The  last  term  is  at  least  of  the  order  of  V2 .  One  can  decompose 
pa{t)  accordingly: 


pa  (t)  =  p1  (i)  +  p11  (0  +  pin  (f )  (103) 

p1  (t)  corresponds  to  the  first  term  of  Equation  102,  p11  (t)  corresponds  to  the  terms  with 
one  kronecker  delta,  and  pIU  (t)  corresponds  to  the  last  term  of  Equation  102.  pl  (t)  can 
be  written  as 
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^ai \a’ i' 


-i(a-a' )wt 


I 

s-0 ,±1 ,±2 


Is  is£’t 

rai;a'i'e 


(104) 


where  the  elements  of  r;*  are  given  in  Ref.  15.  According  to  Equation  104  the  atomic 
population  probability  executes  the  well-known  Rabi  oscillations.  In  p1  (t),  the  effect  of 
the  scattering  on  atomic  population  probabilities  is  simply  to  modify  the  standard  Rabi 
frequency  by  modifying  the  coupling  coefficient  to  the  radiation  field.  The  off-diagonal 
elements  of  p1  (t)  have  spectral  components  at  w,  w  ±  and  w  ±  2£\  The  components  at 
u;  and  u  ±  £'  correspond  to  the  well-known  Mollow  triplet  in  the  resonant  fluorescence. 
The  components  at  w  ±  2£'  are  generated  by  the  collision  process.  Their  amplitudes 
are  proportional  to  T11;0o  and  vanish  if  the  collision  potential  vanishes.  In  the  limit  of 
strong  coupling, 


S  -*oo=> 

Pii.oAt)  — ♦  J  +  (“I)"  \  (/o  -  ^  cos(n«  |  cos  e\t) 

pL, i.  (0  — *  ~  i)  [cos  sin(n*  I  cos  *10 

—  t  ^/o  —  ^  sin^„  sin^2n*  |cos0|t^ 


(105) 


where 


sin  A,  = 


tVinoo  cosfl _ 

<f4[E1  -  E0)7  +  |V11;0o  I2 


(106) 


These  are  as  one  might  have  expected  on  the  basis  of  the  Rabi  oscillations  except  for 
the  second  term  in  the  off-diagonal  element.  pn  (t)  and  pIlI(t )  contain  the  full  spectrum 
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of  pc(t)  :  0,  ±£',  T2£',  w,u;  ±  £',w  ±  2£',2u/,2 u  ±  £' ,  and  2w  ±  2£\  All  of  the  compo¬ 
nents  of  pn  {t)  and  pIIl(t)  vanish  as  the  collision  potential  vanishes.  The  appearance  of 
these  spectral  components  as  a  result  of  the  collision  process  suggests  that,  for  strong 
potentials,  the  Mollow  spectrum,  as  well  as  the  spectrum  of  the  Rabi  oscillations,  is  sig¬ 
nificantly  altered  by  collisions. 

Finally,  the  density  matrix  solutions  are  sensitive  to  the  initial/boundary  conditions.  At 
the  exact  threshold  for  atomic  population  inversion,  both  pa(t)  and  pc(t)  become  diago¬ 
nal  and  time-independent. 

2.3  NONLINEAR  LIGHT  PROPAGATION 

The  preceding  discussion  was  confined  to  a  single  collision  complex  in  a  coherent  field. 
The  following  discussion  concerns  the  problem  of  N-collision  complexes  in  a  coherent 
field,  where  N  is  a  very  large  number.  The  complexes  are  identical  but  distinguishable. 
They  are  initially  independent  of  each  other.  A  quantum  theoretical  treatment  of  how 
a  collective  mode  of  the  complete  system  of  the  coherent  fundamental  model  plus  the 
N-complexes  is  formed  at  the  second  harmonic  is  presented. 

The  standard  semiclassical  theory  of  two-wave-mixing  in  gases  uses  a  physical  picture  in 
which  plane  waves  proceed  from  one  atomic  or  molecular  system  to  another.  To  obtain 
a  measurable  wave-mixing,  the  phases  of  these  waves  must  match  over  a  macroscopic 
distance,  which  requires  the  analysis  of  indices  of  refraction  at  various  frequencies.  The 
following  analysis  precludes  the  use  of  the  semiclassical  theory  for  two  reasons:  (1)  The 
preceding  analysis  used  the  quantum  theory  and  a  representation  which  mixes  the  mat¬ 
ter  and  radiation  states  for  one  complex.  To  be  able  to  use  the  preceding  results  in  a 
straightforward  way,  one  must  remain  within  the  rigorous  picture  of  the  quantum  the¬ 
ory.  (2)  The  single  complex  solutions  were  obtained  for  the  straightline  trajectory  ap¬ 
proximation  in  the  momentum  space  representation  (hence  the  photon  momentum  was 
suppressed).  It  is  necessary  to  provide  a  quantum  picture  of  light  propagation  and  phase 
matching.  This  is  precisely  the  quantum  theoretical  problem  of  finding  a  specific  collec¬ 
tive  excitation  of  a  complicated  many-body  system. 
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The  scattering  states  forming  the  basis  of  the  Hilbert  space  have  precisely  defined  en¬ 
ergies.  Therefore  the  scattering  time  for  these  states  is  indefinitely  large.  These  states 
are  also  in  the  momentum  space  representation.  For  a  scattering  state,  the  spatial  ex¬ 
tents  of  the  particle  trajectories  fill  the  entire  volume  of  the  gas.  Thus,  the  phase  of  a 
state  of  the  complex  plus  the  radiation  field  is  determined  by  the  contributions  of  all 
wavelengths  that  can  be  filled  in  the  entire  volume  and  of  all  points  on  atomic  trajecto¬ 
ries.  This  is  illustrated  in  Fig.  (la).  If  another  pair  of  atoms  is  added  to  the  system,  the 
picture  shown  in  Fig.  (lb)  should  be  superposed  on  Fig.  (la).  For  N  pairs,  there  are  N 
pictures.  The  collective  excitations  are  given  by  what  survives  the  N-fold  superpositions. 
Consider  the  relationship  between  a  macroscopic  property  of  the  N-complex  system  and 
the  corresponding  microscopic  property  of  one  of  the  complexes.  Of  particular  interest  is 
the  transverse  polarization,  which  is  additive. 

Without  the  radiation  field,  a  quantum  state  of  the  N  independent  complexes  can  be 
written  as 


ww».>>=n  w-  (io7) 

p=  i 

where  |/ip)+  designates  the  scattering  state  of  the  pth  complex  and  expfi^  (sp)]  is  an 
arbitrary  phase  factor.  Since  the  complexes  are  independent,  relative  phases  of  their 
states  are  arbitrary.  Here  we  made  the  phase  <f>^  of  the  pth  complex  in  state  |/x)+  de¬ 
pend  on  a  stochastic  variable  (or  a  set  of  variables)  sp .  As  one  goes  from  one  complex 
to  another,  sp  varies  randomly.  Note  that  s  becomes  the  angle  between  e  and  <T10  for  a 
collision  complex  characterized  by  Equations  79-83.  Equation  107  describes  a  multipar¬ 
ticle  state  in  which  each  complex  is  in  a  specific  state  \fi)+  and  the  multiparticle  state 
corresponds  to  a  particular  configuration  An  arbitrary  state  of  the  N-complex 

system  will  be  a  superposition  of  such  states: 


l*>=  £  C(^...nn) (108) 

Mi-Mu 


33 


The  elements  of  the  density  matrix  for  the  pth  complex  (in  the  scattering  state  represen¬ 
tation)  are  determined  from  the  amplitudes  C  and  the  phases  <f> : 


► 


pc'?  = 


C*  (Ml-Mp-lWip+lMn) 


x  1  -pn) 


(109) 


When  a  coherent  radiation  field  jn)  is  introduced,  the  state  of  the  total  system  at  t  =  0 
can  be  written  as 


|*n)t=0  =  •••**») 

°° 


Ml— *»» 

”)«-/*  V'&tlln)®|/xp)+] 


p=l 


n  =  0 


\/n! 


(110) 


where  the  state  of  each  complex  is  associated  with  the  radiation  state  |n). 


Consider  now  the  evolution  of  this  system  with  the  total  Hamiltonian 


v  =  '£hm 


(hi) 


From  Equations  36,  60  and  61,  one  has 


(n  -  m,)/i)+  =  |n)  ®  \p.)A 


(112) 


I")  ®  |(i)+  =  I&0 


(113) 


where  mu  is  defined  in  Equation  71.  Therefore, 
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(114) 


AT 


|^»)t=o  =  £  C(^i-^»)II 

Hi-. .fin  p=l 

(n)n/3 


n  vn!  ^ 


Application  of  exp[— *)/  t/ft]  on  j^n  )t=0  and  the  transformation  back  to  the  original 
Hilbert  space  [using  Equation  61]  yields 


_  JV 

I *n)t=  £  JJ 

i»i  •••»**  p=i 


X  f„AM'(sp)e 


(5 ) 


Am1 


•  £»*(«, )«/*  |n  _|_  m<4,^  ^  |j*'p)+] 


(115) 


Let  A  be  an  operator  which  describes  an  additive  property  of  the  N  complexes: 


*  =  £ap  (116) 

p 

where  Av  represents  the  corresponding  operator  for  the  pth  complex.  The  expectation 
value  of  A  for  the  state  (Equation  115)  is  given  by 


t(^n|X|^n)t  =  £  C‘(n1...pn)C(n'l...fi'n) 


x  £{  n  [e'<^<*»)+<«")e-n£^-£ 


x 


m  M  ,  A  ji 


M‘"n- 


m'  ,Xu' 
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X 


{■ 


Vnln7! 


E 

1  1  ' 


^  g*[^  »x  (*  »o  )  ~  ^  *' x‘  (*»o  M*/*1  ^  ^  &n  + 


m  ii  ,«'  +  m 


MM 


x  ^n-m,1,i)i(spo)^rii(i"('SPii)+  0*po  I^Po  1^^  + 

*  C'-m^i.iV  (5Po  )ln'\'n"'  (Sp0  )  ^  (H^) 

The  first  curly  bracket  reduces  to  a  product  of  Kronecker  deltas.  This  can  be  seen  from 
the  orthogonality  of  the  states  \rpnx)-  First,  from  Equations  36,  59  and  61,  one  has 


(0«a|V»»'a'>  =  5^Cam,»,*V+(m|/*'>+ (»  +  »»#.  K  +  ”V> 

mm' 

n '  ^  iXii^nX' H  ^rni'^H1 


(118) 


Second,  from  the  identity. 


^nn'^M'  =  «n'l  ®  +<M'|)  ’  (|»)  ®  |m)+) 

~  )  1  ^n'-m^x'u'^n- hi^.Am  A'  l^n*) 

A  A  ' 

X 


one  obtains 


^^-m^i.iM^n-pn,,*)!  ~  ^MM'  (120) 

A 
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Thus, 


N 


P*Po 


N 


X  ^-m»1Xn(5p)^»-in|ll,ip'(Sp)  —  PJ 


P#Po 


(121) 


The  expectation  value  (Equation  117)  becomes 


E{E«(o)'-“ 

o=l  '  nil' 


(n) 


Vnln'l 


X  K'-m^.X'ii'  (Sp)^n'X'M"(Sp)+  (Mp  |A>  |/*p  )+  | 


(122) 


The  expression  in  brackets  is  the  trace  of  the  time-dependent  density  matrix  of  the  pth 
complex  with  the  operator  Ap : 


t(*n\A\il>n)t  =  2^[^  +  (m|pc(5p>0Im')+  +  (m'|^p|m)+1 

P=  1  pm' 

=  X>r[P‘(sp,t),ip]  (123) 

P 

t(^n|i(|^n)+  is  a  functional  of  the  configuration  of  the  stochastic  variables  [sj  s2...sN  ]. 
Each  configuration  occurs  with  a  certain  probability.  To  find  the  macroscopic  value  of  A 
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at  t,  an  average  over  all  possible  configurations  must  be  taken.'  For  the  pth  complex,  let 
this  stochastic  average  be  designated  by 

/(<*»,)/(»,)  =  (/).  (124) 

J (dtp)  •  1  =  1  •  (125) 

Thus 


(t{*fi\A\*h)t).  =  fl  / (dsp)t(**\A\**)t 

p=  i  J 

=  N j  (d*)Tr{p'(s,t)A\ 


(126) 


This,  now,  is  the  method  for  obtaining  the  collective  excitations  of  the  entire  system 
directly  from  the  solutions  for  one  complex  plus  the  fundamental  mode.  The  stochas¬ 
tic  average  of  Equation  126  produces  the  phase  matched  solutions  which  can  propagate 
through  the  N-complex  +  the  radiation  fields. 

2.4  SECOND  HARMONIC  GENERATION  FROM  COLLISION  COMPLEXES 
Consider  the  expectation  value  of  the  operator  D: 


0(t)}=  ^2  (127) 

nAn'A' 

where  (*)  giyen  by  Equations  74  and  75.  In  Equation  127,  the  angle  9  between 

the  atomic  dipole  moment  operator  d10  and  the  polarization  vector  e  is  the  stochastic 
variable.  Now  use  Equation  126  and  designate  the  result  by 


Hi)  =  7/"  dHD(t)) 

*  Jo 


(128) 
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P(t)  represents  the  macroscopic  polarization  of  ther  N-complex  system.  The  anomalous 
SHG  arises  from  the  spectral  components  of  P{t)  which  are  at  2u>,2w  ±  £'  and  2 w  ±  £'. 
For  an  intense  coherent  field  and  to  lowest  order  in  V ,  the  amplitude  of  the  component 
at  2 w  is  [Ref.  23]: 


PsH  (0  ~e  e  3,Wt  "rirT^(^m^OO;lo)(^10;ll  —  "^00;0 1  )& (?)  +  C‘C- 


(129) 


where  f  is  given  by  Equation  93  and 


5(f)  m  -f/(l  +  f2)  4-  (3/f)  -  (3 l<?)tan's  .  (130) 

Note  that  in  Ref.  23  a  slightly  different  basis  was  used,  which  introduced  an  extra  y/2 
factor  in  the  coupling  coefficients.  The  results  are  the  same  otherwise. 

Figure  2  shows  S  versus  f.  For  extremely  small  f,  S  is  proportional  to  f.  For  finite  f,  S  in¬ 
creases  at  first  as  f2 ,  then  as  f .  To  the  right  of  the  peak,  S  is  proportional  to  f 1  and 
thus  decreases  as  f  increases,  f2  is  proportional  to  the  fundamental  intensity  Iw .  S2  is 
proportional  to  the  second  harmonic  intensity  /2u) .  Figure  2  suggests  that  in  terms  of 
intensities, 

P  forfl/j  <  10”  2  A 
Pu  fla-lO^A 
Iu  CIr  <  A 

i~  •  nR  >  a 

To  be  able  to  compare  these  results,  for  example,  with  the  Na- vapor  experiments,  several 
complications  need  to  be  addressed.  In  the  atomic  structure  of  Na,  there  are  a  group 
of  states  which  have  approximately  twice  the  energy  of  the  first  excited  states.  This 
fact  has  two  influences.  One  is  that  in  the  calculation  of  +  (0|O|3)+  ,  they  cause  a  near- 
resonant-enhancement.  The  other  is  that  they  can  directly  participate  in  radiative  tran¬ 
sitions,  coupling  both  to  the  second  harmonic  mode,  as  well  as  to  the  fundamental  mode 
(via  two-photon  transitions).  In  other  words,  the  set  of  basis  states  in  Equations  23-26 
must  be  enlarged  to  analyze  the  experiments  correctly. 


(131) 
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The  near-resonant-enhancement  of  +  (0|.D|3}+  may  be  estimated  in  the  following  way. 

The  relevant  group  of  higher  excited  states  may  be  simulated  in  the  present  model  by 
assuming  that  there  is  a  second  excited  a  =  2  state  such  that  E2  «  2 Ex ,  but  the  two 
energies  E2  and  2 Ex  are  sufficiently  different  so  that  there  is  no  need  to  enlarge  the 
set  (Equations  23-26)  within  the  rotating-wave-approximation.  For  the  calculation  of 
+  (0|D|3)+  ,  the  factor 

(^10;l  1  —  Voo-.Ol) 

in  Equation  129  should  be  replaced  by 

^03;11  ^10, 03  +  ^10;30  _  ^00;13 

(2E1  ~  E2 )  _  Ei  Ei  +  E2 

The  comparison  of  the  present  theory  with  the  Na-experiment  is  further  complicated 
by  the  fact  that  the  experimenters  attempted  to  increase  SHG  by  tuning  to  a  double¬ 
photon  resonance  of  the  isolated  atom.  The  results  of  these  attempts  might  have  been 
to  induce  real  population  transfers  to  certain  single  atom  states,  creating  favorable  con¬ 
ditions  for  stimulated  emission  from  these  states,  as  well  as  causing  the  vapor  to  act  as  an 
amplifier  to  the  fundamental  mode.  As  discussed  in  Section  2.2,  nearing  a  population 
inversion  destroys  the  coherence  of  the  collision  complexes.  Without  these  coherences, 
there  is  no  collision-induced  SHG. 

Nevertheless,  with  the  above  caveats,  it  is  possible  to  make  preliminary  comparisons. 

The  best  data  for  this  purpose  is  provided  in  Ref.  6.  The  authors  noted  that  the  second 
harmonic  intensity  varied  as  II  at  first,  then  1% ,  and  at  the  highest  intensities  in  their 
experiment,  as  /„ .  This  is  precisely  the  behavior  described  in  Equation  131.  Clearly,  in 
their  experiments,  they  were  confined  to  an  extremely  small  region  near  the  origin.  In 
fact,  for  the  data  they  label  as  “3s  -  4d”  tuning,  they  did  not  leave  the  region  £exp  < 

0.07,  where  I2u  ~  1^ .  It  follows  that  their  maximum  S  was  about  Ses  «  0.005.  Figure 
2  implies  that  the  maximum  possible  S  is  Sm  «  0.3.  Their  observed  quantum  efficiencies 


(*?.*)  were  on  the  order  of  10" 4  —  10~5.  Therefore,  the  present  calculations  suggest  that 
the  maximum  possible  quantum  efficiency  in  the  Na-vapor  is  on  the  order  of 


*?m  = 


-0.04. 


(132) 


2.5  A  SUITABLE  SYSTEM  FOR  THE  DOUBLING  OF  THE  OXYGEN-IODINE 
FREQUENCY 

The  preceding  discussion  on  Na  points  out  two  important  features  of  the  collision- 
induced  SHG.  The  first  is  that  the  fundamental  frequency  should  be  appropriately  de¬ 
tuned  from  the  atomic  transition  frequency  which  is  to  be  exploited  for  SHG.  This  de¬ 
tuning  depends  on  the  intensity  of  the  fundamental  mode.  The  second  is  that  there 
should  be  higher  excited  states  at  twice  the  fundamental  frequency,  which  contribute 
to  near-resonant-enhancement  of  the  collisionally-induced  coupling  coefficient  to  the  sec¬ 
ond  harmonic  mode.  The  latter  feature  is  not  a  common  property  of  atomic  structures. 

It  is,  however,  much  more  easily  obtained  in  gas  mixtures  in  which  collision  complexes 
of  interest  are  formed  by  two  different  types  of  atoms.  In  particular,  mixtures  of  I  atoms 
with  La  atoms  might  provide  an  efficient  second  harmonic  generator  at  1.3//. 

Although  the  preceding  discussion  was  confined  to  one  type  of  atoms,  certain  qualita¬ 
tive  conclusions  may  be  drawn  from  it  for  an  optimal  system  of  SHG  for  the  1.3//  radi¬ 
ation.  An  analysis  of  the  collision-induced  coupling  coefficient  to  the  second  harmonic 
mode  in  Sections  2.1  and  2.4  shows  that  if  the  collision  complex  is  composed  of  an  A 
atom  and  a  B  atom,  then  there  is  a  term  in  the  coupling  coefficient,  whose  denomina¬ 
tor  is  given  by  (Ef  -  E£  4-  E°  -  Ef),  where  E*'°  is  the  ground  state  energy,  E*  'n 
is  the  energy  of  the  first  excited  state,  and  E*  B  is  the  higher  excited  state.  This  means 
that  for  near  resonant  enhancement  of  the  coupling  coefficient,  the  sum  of  the  transition 
energies  of  the  first  excited  states  of  the  A  and  B  atoms  must  be  near  the  transition 
energy  of  the  second  excited  state  of  the  B  atom.  It  is  possible  to  make  this  denomina¬ 
tor  quite  small  by  choosing  the  I  atom  for  A  and  the  La  atom  for  B.  The  ground  state 
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of  the  La  atom  is  a  *D3/7  state.  It  has  a  *  F£/2  state  which  is  above  the  4  D3/2  state  by 
about  7,000cm-1.  Above  the  state,  it  has  a  bunch  of  states:  4F°/3[l3,600cm-1], 
4D3/2[14,700cm-1],  4£>?/3[14,000cm-1],  4F?/3[15,0G0cm-1  J,  etc.  Here  the  indicated 
energies  are  relative  to  the  ground  state  of  the  La  atom.  For  the  I  acorn,  Ef  —  E£  c* 

7, 600cm" 1 .  Therefore, 


(Ex  -  E0)  +  -  Ef)  as  14,600cm-1 

This  value  is  extremely  close  to  the  transition  energies  of  the  higher  excited  states  of  the 
La  atom  listed  above.  One  would  therefore  obtain  a  strong  resonant  enhancement  of  the 
coupling  coefficient. 

To  be  able  to  provide  quantitative  estimates,  one  needs  to  extend  the  theory  of  the  pre¬ 
vious  sections  to  gaseous  mixtures. 
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3.0  NEW  CRYSTALS  FOR  FREQUENCY  UP-CONVERSION 


3.1.  BAND  STRUCTURE  OF  SUPERLATTICES/MULTIPLE-QUANTUM-WELLS 

In  heterostructure  crystals,  such  as  quantum  wells  and  super  lattices,  nonlinear  suscepti¬ 
bilities  are  increased  several  orders  of  magnitude  over  those  of  the  composite  materials 
in  the  homogeneous  crystal  phase.  They  are  therefore  quite  attractive  for  frequency  up- 
conversion,  and  in  particular  for  harmonic  generation.  As  a  preliminary  to  investigations 
of  their  optical  interactions,  the  principal  investigator  developed  a  new  band  theory  of 
superlattices  and  multiple-quantum-wells  (SL/MOWs),  whose  details  of  the  derivations 
are  given  in  Refs.  24  and  25. 

Since  it  became  possible  to  grow  SL/MQW  crystals  in  the  early  seventies  [Ref.  26], 
there  have  been  several  theoretical  approaches  to  their  band  structure  [Refs.  27-40]. 
These  approaches  are  generally  based  on  some  form  of  either  the  envelope  function  ap¬ 
proximation  [Refs.  33-36]  or  the  tight  binding  approximation  [Refs.  37-40].  The  early 
empirical  results  concerning  absorption,  luminescence,  carrier  transport,  etc.  in  super- 
lattices  could  be  explained,  at  least  qualitatively,  by  means  of  simple  quantum  wells  and 
Kronig-Penney  models  [Refs.  27-32].  Later  experiments  revealed  the  appearance  of  cer¬ 
tain  forbidden  transitions  [Refs.  41-42].  To  explain  these,  as  well  as  subband  energy  dis¬ 
persions,  one  needs  to  go  to  more  sophisticated  theories  in  which  band  mixing  occurs 
[Refs.  35,  36,  43]. 

In  Bastard’s  envelope  function  approach  [Refs.  33-36],  one  generalizes  Kane’s  theory 
[Ref.  44]  for  the  IH-V  semiconductor  compounds  to  superlattices,  and  writes  a  Kane 
Hamiltonian  for  each  type  of  material  composing  the  superlattice.  The  solutions  of  these 
effective  Hamiltonians  are  then  joined  across  layer  interfaces  by  means  of  boundary  con¬ 
ditions  which  require  continuity  of  wave  functions  and  conservation  of  electron  current 
densities.  These  calculations  ultimately  depend  on  the  effective  mass  concept.  The  ef¬ 
fective  mass  of  the  electron  is  assumed  to  be  position  dependent  along  the  superlattice 
axis  (that  is  to  say,  along  the  direction  which  is  vertical  to  the  planes  of  the  superlattice 
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layers) .  The  space  dependence  of  the  effective  mass  is  not  calculated,  but  is  postulated 
to  have  a  certain  form;  then  the  results  are  compared  with  empirical  observations  to  the 
extent  that  this  is  possible. 

Tight  binding  calculations  of  superlattices  can,  in  principle,  be  more  accurate  than  cal¬ 
culations  based  on  the  envelope  function  approximation,  but  they  involve  extensive  and 
complex  numerical  computations.  In  a  tight  binding  calculation,  one  uses  superpositions 
of  large  numbers  of  electronic  orbitals  to  solve  for  the  eigenvalues  and  eigenfunctions  of 
a  tight  binding  Hamiltonian  [Ref.  45].  Band  mixing  can  be  directly  incorporated  into 
these  computations  and  optical  spectra  of  superlattices  can  be  reproduced  reasonably 
accurately  [Ref.  43]. 

This  report  presents  a  new  one-electron  theory  of  superlattices.  The  homogeneous  crys¬ 
tal  band  structures  of  the  composite  materials  are  assumed  to  be  completely  known.  A 
superlattice  Hamiltonian  in  the  momentum  Bloch  representation  is  obtained,  which  ex¬ 
presses  all  of  the  effects  of  symmetry  breaking  in  terms  of  the  quantities  and  parame¬ 
ters  related  to  the  homogeneous  crystals  of  the  composite  materials  and  the  superlat¬ 
tice  geometry.  Also,  series  expressions  for  the  superlattice  subband  energies  and  wave 
functions,  which  can  be  evaluated  systematically  with  increasing  accuracy  are  presented. 
The  virtue  of  the  new  formulation,  however,  lies  in  the  fact  that  it  provides  much  greater 
conceptual  clarity  and  pliability  for  analytic  manipulation  than  the  previous  theories. 

Consider  two  types  of  materials,  A  and  B,  whose  homogeneous  lattice  structures  match: 


RA  —  Rg  —  R 


(lattice  translation  vectors) 


(133) 


Ga  =  GB  =  G  (reciprocal  lattice  vectors) 


(134) 


Assume  that  the  Bloch  functions  of  the  homogeneous  crystals  of  these  materials  are 
completely  known: 


B*b  =  2^  +  Vab(£> 


(135) 
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HA,B\nk)A.B  =  E*'B  {k)\nk)A'B-,(x\nk)A'B  =  ip*j.D(x) 


(136) 


Assume  that  the  reciprocal  lattice  is  a  simple  cubic  lattice  with  one  of  the  principal  di¬ 
rections  along  the  z-axis.  This  is  a  minor  simplifying  assumption;  the  results  are  easily 
generalized  for  arbitrary  lattice  orientations.  The  reciprocal  lattice  vectors  Gz  can  be 
written 


G,  =  —M,  M  =  0,±1,±2, .  (137) 

a 

a  is  the  fundamental  lattice  constant  of  the  homogeneous  crystals  of  A  and  B.  Define 

l0  =  —  =  a  positive  integer.  (138) 

a 

Consider  now  the  geometry  shown  in  Fig.  3.  The  layers  of  material  A  of  width  LA  are 
sandwiched  in  between  the  layers  of  material  B  of  width  LB  .  Assume  initially  that  a 
finite  number,  M,  of  A  layers  are  embedded  in  B.  In  the  momentum  Bloch  representa¬ 
tion,  let  U  become  infinite,  making  the  composite  structure  a  superlattice.  Keeping  M 
initially  infinite  avoids  some  minor  mathematical  problems,  as  discussed  in  Ref.  24.  Let 
the  z-axis  be  vertical  to  the  layers  (the  superlattice  axis).  Pick  the  origin  of  the  coordi¬ 
nate  system  as  symmetrically  as  possible  for  the  convenience  of  evaluating  certain  sums 
later.  If  U  is  an  odd  integer,  the  origin  is  in  the  middle  of  an  A  layer  as  shown  in  Fig. 
3a.  If  N  is  even,  the  origin  is  in  the  middle  of  a  B  layer  as  shown  in  Fig.  3b.  In  the  su¬ 
perlattice  limit,  whether  one  starts  with  odd  or  even  M  does  not  matter;  the  results  are 
independent  of  the  initial  configuration.  For  the  following  choose  the  initial  configura¬ 
tion  to  be  the  one  shown  in  Fig.  3a  with  odd  XI.  The  layers  of  A  and  B  are  assumed  to 
have  perfectly  sharp  plane  interfaces. 
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(a) 


Origin  of  the 


coordinates  for  odd  (a)  and 


The  one-electron  Hamiltonian  of  the  composite  crystal  of  Fig.  3a  can  be  written  as 


H  =  +  V0  (*)  +  |V„  (i)  -  VB  (i)|  E'.;  I  Ji'I  „„ 

x  ffl(2-sL+ii^)-«(2-sL-iiA  =  ~+VB(z)  +  Htw(x)  (139) 
2  2  2m 

where  0’s  are  the  step  functions,  and  L  is  the  superlattice  period  along  the  superlattice 
axis: 

L  =  La  ■+■  Lb  (140) 

Consider  now  the  matrix  elements  of  this  Hamiltonian  in  the  representation  defined  by 
the  Bloch  functions  of  the  B  type  crystal.  In  the  superlattice  limit,  namely  when 

U  — ►  oo 


(141) 


(142) 


a  (nk[H,„  \nk')B  =  aa.f(k,  k‘  +  —z) 

x  W(G'  -  6)4° '(k  -  G)4°Ak-G  ~~i+  ~i)  [143) 

where  <t>*  (A:)  is  the  momentum  Bloch  function  of  the  crystal  type  B  and  W  (G)  is  related 
to  band  offsets: 


a,-„[i^(l-loM)1 


W(G)  =  V|  -  V| 

=  E„  (0)«  (0)^'  (-G)  -  E ;  (0)tf  (0)*f  (— G)]  (144) 

Next,  consider  a  superlattice  eigenstate  \tp).  Since  the  set  of  B  type  Bloch  states  is  com¬ 
plete, 


W  =  S„Jan(i)|nf)J,  (145) 

Note  that  the  sum  over  k  is  confined  to  the  Brillouin  Zone  of  the  homogeneous  crystals. 
The  time-independent  Schrodinger  equation  for  |t/>)  is 

H\1>)  =  Xn,£,an.{k'){HB  +  HAW)\n'k')B  =  c|0)  (146) 

Define 

k  =  k  +  ^-z  (147) 

Li 

such  that 

kx  =  kx 


=  k„ 


7T  7T 

~l<k-<l 


1=0,±1,±2, . ±(1„-1) 


(148) 
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.  _  27 r/ 

<x„i(/c)  =  an(/c  +  —  2) 


(149) 


«;.(*)•=  e;  (« *) 


(150) 


A  -!0M]r  TirtsZi 

v.„„.  (g) = ,l*  sr—  ‘  y_ -  foM]  Sg°’lv(G  - 


x  <£ f  («  —  G  + 


-  27rM . 

G  + - 2  + 

a 


2tt/'  \ 

~ ' V 


Equation  145  becomes 


(151) 


[c  —  E%t  (*)]  ani  (A)  —  £„<!'  Vnl-.n'l'  (*)  ~  0 


(152) 


««(3  = 


anJ0-l(«) 

Cnlo-2  (^) 
an-<0  +  1  (*) 


(153) 


(*)  = 


(154) 


l^nr,'  (*)]„<  =  Vnl.n'l'{> «) 


(155) 


Equation  152  can  then  be  written  in  terms  of  the  2/0  -  1  component  spinors  an  and 
(2/0  -  1)  x  (2 10  -  1)  matrices  Vnn-  as 


£„•  [tn  +  Kn' (*)]  <*»'(*)  =  «<*»(*) 


(156) 
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Equations  152  and  156  show  that  a  given  band  of  the  homogeneous  crystal  splits  into 
2/0  —  1  components  in  the  superlattice.  The  components  are  the  superlattice  subbands. 
There  is  a  correspondence  between  them  and  the  bound  states  of  quantum  wells.  The 
splitting  arises  entirely  from  the  symmetry  breaking  along  the  i-axis.  The  original  Bril- 
louin  Zone  of  the  homogeneous  crystal  is  partitioned  into  2/0  —  1  sections.  The  central 
piece  for  —  <  /c,  <  f  becomes  the  new  Brillouin  Zone  for  the  superlattice.  The  sub¬ 

bands  of  a  band  n  arise  from  the  folding  of  the  original  En  onto  the  superlattice  Bril¬ 
louin  Zone.  Furthermore,  the  Hamiltonian  which  determines  the  subband  energies  and 
subband  states  (in  the  representation  of  (|n£) B  })  becomes  a  supermatrix  given  by 


+  v„,„ 


'(*) 


(157) 


Note  that  km  =  ±ir/L  represent  identic^  1  points  in  the  superlattice  Brillouin  Zone,  since 
they  are  related  by  the  momentum  2 ttz/L.  Only  one  value  is  needed  in  the  superlattice 
Brillouin  Zone,  for  example  kx  —  +nlL  . 

Simpler  expressions  than  Equation  151  can  be  obtained  for  Vnl.n>r  with  the  aid  of  the 
interband  position  operator.  They  are  given  in  Ref.  24. 

The  energy  eigenvalues  and  eigenstates  of  Equation  152  may  be  obtained  using 
Lowdin’s  theorem  [Ref.  46].  Define  the  matrix  function  (E)  in  terms  of  the  se¬ 

ries 


V  V' 

2(  jr»\  _  ySL  I  rv  nl\n\li 

filin' l*  [A)  “  *nl,n'l'  +  ^nilx 


"1" 


V' 


r»|  njljin'l' 

(E  -  X?ih)(E  - 


+  •  *  . 


(158) 


where 


yD  _  ySL 

nnl  —  *nl;nl 


(159) 
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(160) 

The  primes  over  the  summation  signs  mean  that  the  sums  exclude  the  ni- terms.  An 
eigenvalue  of  Equation  152  is  given  by  the  solution  of  the  equation 


Vnt;n'l'  ~  (l  —  &nn'&ll,)Vnl-,n'l 


en,Z)  =  ii.i(£„i(ic)) 


(161) 


such  that 


lim  eni (k)  =  E^{k) 
w  ->  o 


(162) 


The  nondegenerate  eigenstate  corresponding  to  enl  is  given  by 


V'm*  (*)  =  r/n,  < 


Zn't‘;nl{en,)rl>*  ^(X) 
,  yi  _ n  K+  l  * 


(163) 


where 
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Reference  24  describes  the  procedure  in  case  of  degeneracy. 


To  explore  the  implications  of  the  preceding  results,  let  us  consider  a  superlattice  crystal 
composed  of  two-band  materials  (one  conduction  and  one  valence  band).  Assume  that 
these  bands  are  parabolic  in  the  homogeneous  crystal  phase  and  that  there  is  no  degen¬ 
eracy.  With  further  approximations  discussed  in  Ref.  24,  the  superlattice  conduction 
subbands  are  given  by 


eeo{Z)  =  E?0{K)  +  pLAW{  0)+E« 


[£?„(*)-£=(*)]  (4-1)' 

jin^ 
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for  the  lowest  subband  l  =  0,  and 


1 


for  finite  l.  Here  £  is  a  band  mixing  parameter  [see  Equation  80  of  Ref.  24]. 

Consider  given  by  Equation  166.  The  first  two  terms  define  the  position  of  the  en¬ 
ergy  band,  taking  into  account  the  band  offsets  between  the  two  materials.  The  third 
term  gives  mainly  the  dispersion  in  the  transverse  directions.  The  fourth  term  comes 
from  the  folding  of  the  Brillouin  zone  and  describes  the  partial  localization  of  the  elec¬ 
tron.  For  LA  =  Lb  =  Lj2,  this  term  becomes 


.  r  ft2*2  2n2h2l2 
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This  is  identical  to  the  standard  expression  for  a  rectangular  well  of  width  LA  (for  an 
infinitely  deep  well). 

The  above  comparison  with  the  quantum  well  model  depends  on  the  fact  that  the  homo¬ 
geneous  crystal  bands  Me  parabolic.  If  they  are  not  parabolic,  it  is  much  harder  to  make 
a  comparison.  For  some  narrow  gap  semiconductors,  conduction  and  valence  bands  can 
be  approximated  with  linear  dispersion  relations.  For  example,  let 


E?  (k)  =  Ef  (0)  +  \a  •  k\ 

Conduction  subband  energies  for  finite  l  become 


(168) 
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where  W,Ct  and  S(  are  defined  in  Ref.  24.  The  dependence  on  |/|  here  is  quite  different 
than  in  Equation  167. 

The  last  term  of  Equation  166,  which  is  a  square  root,  makes  the  subbands  nonparabolic 
in  the  z-direction.  In  the  x,  y-directions  they  are  still  parabolic.  Because  of  nonparabol- 
icity,  one  must  be  careful  in  assigning  an  effective  mass  for  a  subband.  The  effective 
►  mass  can  be  quite  different  for  «,  =  0  and  kx  =  finite  for  nonparabolic  bands.  For  /  ^  0, 

the  effective  masses  as  km  — ♦  0  axe  given  by 
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1± 


4  n2h2l2 
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For  l  =  0: 
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Note  that 
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mc0 

mBc 


>  1, 
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(172) 


For  large  /,  the  difference  of  m^LJmBc  from  unity  becomes  more  pronounced.  These 
theoretical  predictions  for  the  subband  effective  masses  should  exhibit  themselves  in  lin¬ 
ear  and  nonlinear  susceptibilities  through  the  density  of  states. 

Band  mixing  occurs  throughout  the  formulas.  In  the  subband  energies,  band  mixing  is 
associated  with  the  parameter  £.  If  £  is  zero  or  some  multiple  of  i r,  then  band  mixing 
does  not  occur.  In  the  subband  wave  functions  given  by  Equation  163,  the  last  set  of 
terms  under  the  summation  sign  represents  the  band  mixing  terms.  These  terms  are  re¬ 
sponsible  for  the  modification  of  the  selection  rules  for  radiative  transitions. 
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When  the  electronic  spin  is  taken  into  account,  the  Hamiltonian  (Equation  142)  becomes 
much  more  complicated.  The  spin-orbit  coupling  is  discussed  in  Ref.  25.  Here  the  dis¬ 
cussion  is  restricted  to  a  few  comments. 

The  spin-orbit  coupling  contributes  two  terms  to  the  overall  Hamiltonian:  Xs  1  [see 
Equation  19d  of  Ref.  25]  and  Xs2  [Equation  19e  of  Ref.  25],  Xs  1  is  related  to  the  pres¬ 
ence  of  the  interfaces.  When  an  electron  crosses  an  interface,  it  experiences  a  pulse  of 
magnetic  field  in  its  rest  frame.  Xs  1  describes  the  coupling  of  the  electron  to  this  field. 
Xs 2  arises  from  the  zone-folding  of  the  difference  of  the  spin-orbit  interactions  of  the  ho¬ 
mogeneous  crystals.  The  present  theory  suggests  that  the  off-diagonal  matrix  elements 
of  Xs  1  is  an  order  of  magnitude  larger  than  the  matrix  elements  of  Xs2.  The  main  in¬ 
fluence  of  Xs  1  is  on  the  selection  rules  for  the  radiative  transitions.  These  selection  rules 
are  modified.  The  energy  shifts  of  the  subbands,  however,  come  mainly  from  Xs2 . 
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